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Abstract 

Despite there are numerous theoretical studies of stochastic differential equations with a symmetric a- 
stable Levy noise, very few regularity results exist in the case of 0 < a < 1. In this paper, we study the 
fractional Fokker-Planck equation with Ornstein-Uhlenbeck drift, and prove that there exists a unique 
solution, which is C°° in space for t > 0 when a G (0, 2]. 
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1. Introduction 


In this paper, we are interested in the regularity of solutions of the generalized Fokker-Plank equation 
for a; G and t > 0 

dtU=I[u]-S7 ■ {f{x)u), m(0,x) = uo(a;), (1.1) 

where X is an operator defined in (12.111 . f{x) : —>■ in the drift term is a deterministic driving force, 
and uq in is nonnegative. We assume without loss of generality that UQ{x)dx = 1. When 

X is the classical Laplacian, Eq (HH) is the classical Fokker-Plank equation. If /(■) fulfills the usual 
assumptions of regularity needed in the Laplacian perturbation theory, it is well-known that the solution 
is smooth due to the regularizing effect of the parabolic operator. If /(•) is less regular, there is a series 
of papers by Lions, et al. devoted to vector field /(•) having bounded divergence and some Sobolev type 
regularity, and it is shown that the regularity in x is sharp for /(x) being in Sobolev spaces instead 
of being Lipschitz or continuously bounded (see la Em). 

The Laplacian is the infinitesimal generator of the Brownian motion, hence the classical Fokker-Planck 
equation corresponds to a stochastic system with Gaussian noises. However, the noise is often non- 
Gaussian, and in this paper, we focus on the symmetric a-stable Levy motion, for which the infinitesimal 
generator is the fractional Laplacian, denoted by with a G (0, 2). Eq (11.111 with X = is called 
the fractional Fokker-Plank equation (see Section 2 for details) 


There are many results (see i, a-m, 


and references therein) on gradient 
perturbation of the fractional Laplacian operator in the frame work of regular vector field, i.e. the well 
posedness and regularity of the solution to the equation ut = A“/^m — / • Vit with uq G and 

/ satisfying some regularity assumptions. Depending on the value of a, there are three case^ First, 
a G (0,1) is known as the supercritical case because A“/^u is of lower order than / • Vu. In 
shown that the solution becomes differentiable with Holder continuous derivatives, if / G 
any /? G (0, a). With a slightly weaker regularity assumption on /, it is shown in 


25 


continuous if / G “in space. Second, in the critical case of a = 1, Silvestre in 


|24| . it is 

that u is Holder 
shows that u 
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becomes Holder continuous for positive time if / S while the authors in Q and [1^ achieved 

the same result by different approaches when / belongs to the BMO class and div/ = 0. Third, when 
a € (1, 2), the term is dominant, which is called the subcritical case, and better regularities have 

been proved under more general conditions. Actually, in the one-dimensional case {d = 1), Droniou, et al. 
in [1^ show that the unique solution exists and u € C°°((0, oo), K) for / ■ Vu replaced by dx{f{u)). Their 
result is generalized to multidimensional cases in [l^ with the term H{t,x,u,Vu) instead of dx{f{u)). 
Note that f{u) and H{t,x,u,Vu) must satisfy certain conditions in order for to be the dominate 

term. The estimate of the fractional heat kernel perturbed by gradient operator is obtained in Q if / 
belongs to a Kato class. See also @,[li for two-sided fractional heat kernel estimates. 

However, for the vector field / lying in some Sobolev spaces instead of being uniformly bounded in 
X, the drift term could be stronger than the diffusion term T[u\ even in the case of I = A. Hence it 
would not be possible to prove any regularity for u in general and there are relatively fewer papers. For 
the fractional Fokker-Plank equation, the existence of weak L^-solution was recently obtained in (^ . 
if /(•) € where q is the conjugate exponent of p. The solution is also unique with further 

assumptions on /(■). However, there is no regularity result so far. 

In the rest of the paper, we focus on the fractional Fokker-Plank equation with the Ornstein-Uhlenbeck 
drift, i.e., Eq (11.11) with X = A“/^ and f{x) = —x. 


dtu = A“/^m — V • (— xm), u(0, x) = uo(a^)- 


( 1 . 2 ) 


In this case, the unique weak L^’-solution exists by applying the results in 2^. However, we show that 
the unique solution to Eq (11.21) is smooth in this work. More precisely, the solution of dOJ can be 
represented as 

(1.3) 


i(i,x) := (p*uo)(x) = / p{t,x - y)uo{y)dy, 


where p{t, x, y) = p{t, y — x) is the fundamental solution of (II.2|) (or p is called the heat kernel of the 
fractional Laplacian perturbed by the Ornstein and Uhlenbeck drift). One can easily check that 


p{t,x,y) = e‘^^p 


- 1 


Q 


e*x,?/ 


(1.4) 


where p{t,x,y) = p{t,y — x) is fundamental solution to the fractional heat equation, i.e., dtU = A°‘^^u 
(or p is called the fractional heat kernel). It is well-known that p(i, •) is for alH > 0 when a € (0, 2) 
(see a)- Due to (jl.4p . the solution of has the same regularity as that the fractional heat equation, 
as shown in the following theorem. 

Theorem 1.1. Let uq S L°°(]R^) and be continuous almost everywhere. The funetion u defined in 
don is the solution to dOl). Furthermore, for T > 0 or T = oo, 

u{t,-) €C°°{W^), vte(o,r]. (1.5) 


We remark that from the stochastic view point, the o-stable Ornstein-Uhlenbeck process was first 
studied in the case of a = 1 in 17|, and the sharp estimates of the mean first exist time from the ball 
was given for 0 < a < 2 in , and it is shown in 211 that the Harnack inequality holds if 1 < a < 2 or 
a < 1 = d. Erom the pseudo differential equations point of view, see [isj for the exponential decay in time 
of the solution to the steady state of the associated stationary equation, and 0 for the applications of the 
Levy logarithmic Sobolev inequality to the study of the regularity of the solutions of the fractional heat 
equation and the asymptotic behavior of the Levy-Ornstein-Uhlenbeck process. Eor numerical results, 
see 141. Finally, we refer to 13, 27, 2^, 3^ for applications of the fractional Ornstein-Uhlenbeck process 
in physical sciences. 

We point out that our result in Theorem 1 1.1 1 is a little surprising, especially for the supercritical case 
(a S (0,1)), because it implies that the fractional Laplacian operator A“/^ has the smoothing effect 
like the Laplacian in the Ornstein-Uhlenbeck drift case. We wish to expand the result to include more 
general drift terms in further studies. The paper is organized as follows: in Section 2, we introduce some 
notations and definitions. And the proof of Theorem ll.il will be given in Section 3. 
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2. Preliminaries 


Let us briefly recall basic definitions of the isotropic a-stable Levy operator and the potential theory 
of the fractional Laplacian (see[I|, 0, H, I^L 

The Levy-Khinchine formula implies that there exists a symmetric semi-definite d x d constant matrix 
a = (cTij), a constant vector b = {bi) S and a nonnegative measure v on such that the Levy 
operator writes: 

I[u](a:) = (ctVm; Vm)-I-( 6; Vu)-I- / {u{x + z) — u{x)— Vu{x) • (2.1) 

Js.<‘ 

where (■; ■) is the inner product. The so-called Levy measure v can be singular at the origin and satisfies 

^({0}) = 0 and J min(l, \zfMdz) < oo and the matrix a characterizes the diffusion (or Gaussian) part 

of the operator (it can be null if cr = 0), while b and v characterizes the drift part and the pure jump part, 
respectively. The characteristic exponent of the Levy process is denoted by '0(C) = (o’C; C) + (b', C) + «(C)i 

where a(C) = / [1 — cos(C ■ z)]i/{dz) is associated with the pure jump part of the underlying process. 
Js.'^ 

The symmetric a-stable Levy process refers to the case 6 = 0, cr = 0, and z/(dz) = Ad,c 


dz 


\d+a 


Ad,a is so chosen that a(C) = ICT for C € (see Q and references therein). 
For a G (0, 2) and 0 G the fractional Laplacian is 


A 2 0(x) := Ad,a, lim / 


0(x + z) — 0(a:) 


\d+a 


dz. 


where 


( 2 . 2 ) 


Note that if a = 2 the Levy operator has no jump part, and ![•] is the Laplacian operator in this case. 
We let p{t, •) be the smooth real-valued function on with Fourier transform 


[ p{t,z)A^'^dz = e t > 0, ^G. 

dRd 


(2.3) 


It follows from (IQ) that p is the heat kernel of the fractional Laplacian and it is well-known that for 
t > 0, x,y G X ^ y, the sharp estimate of p{t, x) is 


p(t, X, y) pz min 


_ y\d+OL ’ 


.-d/a 


(2.4) 


(see @,0,0,0), where the notation pz means that either ratio of the sides is bounded by C £ (0, oo), 
and C does not depend on the variables shown, here t and x. The estimate of the first order derivative of 
p{t,x) is was derived in @, Lemma 5], and we extend it to get the estimate of the m-th order derivative 
of p{t, x) for any m > 0 as follows 


.= L^J 


dTKtA)- ^ Cn\x-yr ^”min|- 

n=0 ^ 


d+2(m-n) 


_ y|tZ+a+2(m—n) ’ 


(2.5) 


where means the largest integer that is less than y- 


3. Proof of Main Result 

Proof. To prove Theorem ll.il we first prove the following result: for uq G L° 
almost everywhere. The function u{t,x) defined as 


u{t,x)= / p{t,y - x)uo{y)dy 


and uo(x) be continuous 
(3.1) 
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is the solution to Ut = A“/^u with u(0, x) = uo{x), and for T > 0 or T = oo, we have 

u{t,-) vte(o,r]. (3.2) 

Firstly, we show that u(t, •) G C'°°(R'*) for all t € (0,T]. Since ug € and p(t,-) G for 

t > 0, p(t, •) * uo is well-defined. By the Young’s inequality for the convolution in (13.IL we have 

V(t,a;) G (0,oo) X ||p(t, •) * wo(a;)||Li(Rd) < ||uo||Loo(Rd). (3.3) 

By the smoothness of p{t, x, y), to show ()3.2I1 . it is sufficient to show that 

|V™p(t ,X — y)\ G L^{R'^), Vm > 0,to < t < T, and x,y G R'^. 

Indeed, by (12.511 . we have 


n=[^J 

n=l^\ 

= I \^-y\ 


d+2{m — n) 


m— 2 n^— 


d + 2{m-n) 


d-\-2(m — 


1=0 JB{x,r) 

"=1t1 

< E *0 

n—O 

< oo. 


n=l^i 

dy+ 


dy 


t 


”=1tJ 1 


:dy 


' B{x,r) 


\x-yr-^^dy+ T 


n=0 


B{,x,rY |x-2/|‘^+“+' 


:dy 


Sine Mo G and p{t, •) G C'°“(R‘’*), the theorem of continuity under the integral sign gives ()3.2F 

The uniqueness of the solution u{t, x) defined in (j3.1ll follows from the Fourier Transform. And it is clear 
that for x'^ at which uo(x) is continuous, we have 

lim u(t, x) = uo(x^). (3.4) 

t—^0,x—>^xo 


Indeed, fix a;° G R'^, e > 0. Suppose mo(x) is continuous at xg- Choose d > 0 such that 

\ug{y) — m(x°)| < e if ly — < d, and y G R’’*. 

Then if \x — x°| < |, we have 


i(t,x) - ito(a;°)| < / p{t,x-y)\ug{y)-uo{x°)\dy 

p{t,x- y)\ug{y) - uo{x°)\dy + 
t 


JB{x^,5) 

<e -|- 2||uo||l= 


B{x°,SY 


p{t,x- y)\uo{y) - uo{x°)\dy 


Ib{xO,sy “ 2 / 1 '^'^“ 


dy. 


In the second term, we have \y — x^\ >5. Hence 


(5 1 

\y-x°\<\y-x\ + \x° -x\<\y-x\ + -<\y-x\ + -\y-x^\. 


Thus \y — x\ > ^\y — x'^\. Consequently, the second term in the last line tends to 0 as t ^ 0+. Hence, if 
\x — x°\ < |, and t > 0 is small enough, \u{x,t) — mo(x°)| < 2e. 

Second, set i = and x = e*x. By (II.4|) . we have 

V™p(t, X, y) = e‘'‘e™*V^p(f, i, y), (3.5) 
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Hence by (12.51) . we have 


"=L^J 

|VXi,x-2/)|d2/<e(‘'+-)‘ ^ 

n=0 

n=l^\ 


R’i 


- y\ 


t 


t~ 


d+2(m-n) 


_ g(d+m)t 


E / 

„=0 


2=LtJ 


< ^ f- 


ti4'2(m —n) 


^ y|(i+a+2(m—n) 

n—0 
n=l^\ 

_y\m-2UdyJ^ I 


dy 


t 


\x-y\d+‘-+- 

f 1 


n=0 


' B{x^r) 


=0 'IB{x,rY 


-y\ 


d+a+m 


< OO 


for any 0 < to < t < T. Hence, the theorem follows. 


□ 
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